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Abstract 

In this paper we study in a large-S expansion effects of substituting spins by 
non-magnetic impurities in two- and three- dimensional Heisenberg antiferro- 
magnets in a weak magnetic field. In particular, we demonstrate a novel mech- 
anism where magnetic moments are induced around non-magnetic impurities 
when magnetic field is present. As a result, Curie- type behaviour in magnetic 
susceptibility can be observed well below the Neel temperature, in agreement 
with what is being observed in La2Cu\- x Zn x O '4 and Sr{Cui- x Zn x )20^ com- 
pounds. 

PACS Numbers: 75.10.Jm, 75.10.-b, 75.40.Mg 



1 



Typeset using REVTpjX 



Recently, there has been increasing interests in the effects of substituting Cu by non- 
magnetic ions (Zn) in high-T c cuprates, where interesting effects were discovered for com- 
pounds in the underdoped regime |I| . In particular, magnetic susceptibility measurement |2j 
and other experiments 0] seem to indicate that local magnetic moments are generated as Cu 
ions were replaced by Zn, both in the underdoped high-T c Yttrium and Bismuth compounds 
||||, and for undoped Lanthanum compound Q. Based on RVB theories of the t — J model, 
Nagaosa and Ng || have constructed explanations on how local magnetic moments can be 
generated in these compounds in the disordered phase of two-dimensional antiferromagnets. 
Their theory is believed to be applicable to the underdoped regime of high-T c compounds. 
However the theory is not applicable to the La compound where Curie-type behaviour is 
observed in magnetic susceptibility at temperatures well below the Neel temperature where 
long-range antiferromagnetic order exists. More recently, similar experiment on the two-leg 
ladder compound Sr(Cui- x Zn x )20z also indicates formation of local magnetic moments 
in the spin-gap phase as ions are replaced by Zn ions ||. The compound was found 
to order antiferromagnetically at low temperature at x ~ 0.07 |7J. Surprisingly, Curie be- 
haviour in magnetic susceptibility remains even below the Neel-ordering temperature [|7J, 
indicating that formation of magnetic moments around non-magnetic impurities is a general 
phenomenon in Neel state of quantum antiferromagnets. 

In the presence of long-range antiferromagnetic order, the low temperature properties of 
clean antiferromagnets can be understood at least qualitatively using a semi-classical theory 
(1/S expansion). In this paper we shall generalize this approach to include the effects of 
non-magnetic impurities in the presence of a uniform magnetic field. Notice that in the 
absence of magnetic field, a similar study has been performed by Bulut et. al. || where 
no 'free' magnetic moments were found to be induced around non-magnetic impurities. We 
shall show that local magnetic moments can be induced by non-magnetic impurities once 
magnetic field is applied on the system through a novel mechanism. Notice that strictly 
speaking, thermal fluctuations destroy long-range magnetic order at two dimension at any 
finite temperature. Thus our analysis at two dimension can only be applied to layered 
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systems like Lanthanum cuprates where the magnetic coupling between different layers is 
much weaker than intra-layer coupling. 

To begin with, we first consider the zero temperature problem of classical spins inter- 
acting antiferromagnetically under a uniform magnetic field B, with a single non-magnetic 
impurity replacing spin at site i . We shall consider the Neel state in the absence of mag- 
netic field to be ordered in the ^-direction, with the uniform magnetic field B applied in 
the +x-direction. The classical ground state in the absence of magnetic field has spins all 
pointing in +z direction for spins in A-sublattice and spins all pointing in —z direction for 
spins in i?-sublattice. In the presence of magnetic field, the spins tilt to the +a;-direction to 
minimize the magnetic energy. Let 9f be the angle tilted away from the z-axis for spin on 
site i on A-sublattice and 9 B be the corresponding angle for spin on site j on 5-sublattice. 
It is easy to show that the classical energy E c i is given in the limit when 9 A ^ are small 
(weak 5-field limit) by 

E d /S 2 = - £ (i- ^ + e f )2 )- B "£0f-B"£of, (1) 

<i^io,j> iy^io j 

where < i ^ i , j > are nearest neighbor sites in the square (cubic) lattice excluding contri- 
butions from the non-magnetic impurity at site Zq and B' = gfisB/S. Notice that the only 
effect of non-magnetic impurity is to remove the spin at site i$ in our theory. We have also 
set the spin-coupling J = 1 in Eq. (p. Notice that the energy expression for E c i is valid 
only to order 0(B 2 ). 

Next we consider the continuum limit of the energy expression (|1|). Introducing sym- 
metric and antisymmetric angle variables 

9 s (x) = S(9 A (x) + 9 B (x)), 

9 a (x) = ^(9 A (x)-9 B (x)), (2) 

we obtain after some straightforward algebra 

d d x 



E, 



< 



(d)9 s (x) 2 - B'S9 s (x) + (V9 a (x)) 2 + 9 -^- (9 s (x)d .V9 a (x) - 9 a (x)a .V9 s (x)) 



+E (3) 



where Eq = —NS 2 d/2, N = number of sites in the system, a a is lattice spacing and d D is a 
radial vector pointing away from i . v a is the volume enclosed by the sphere (circle in 2D) 
with radius |. For an impurity located at B— sublattice, d a — > — d a . 

The various terms appearing in Eq. (|3]) can be understood rather easily. First of all, in 
the limit a — > 0, X{ — > Xj in Eq. (JI|) and the only contribution to E c i would be the first 
two 9 S terms in Eq. (§). It is also clear that in the limit a a — > 0, terms proportional to 9 a 
or 9 2 do not appear in Ed, and the only contribution from 9 a can appear as (V6* a ) 2 only. 
Appearance of non- magnetic impurity at site iq breaks the symmetry between A— and B— 
sublattices and introduces local coupling between the 9 S and 8 a fields. Notice that additional 
coupling between 9 a and 9 S fields will appear if we take into account higher order terms in 
B in our energy expression E c \- However, in the absence of impurities, these terms do not 
break the symmetry between A— and B— sublattices or the symmetry of interchanging 9 A 
and 9 B fields. 

Minimizing E c i with respect to the 9 S field we obtain 

ds{S) = Yd { B ' S ~ ^v^«o- v ^( f )) • ( 4a ) 

Putting 9 S back into E c i and minimizing with respect to 9 a field, we obtain for \x — xq\ > a , 

V 2 9 a (x) = ^^5 d (\x-xl\-a ). (4b) 
dv a 

Notice that in the presence of magnetic field, the non-magnetic impurity acts as a source 
term for the 9 a field of strength 2B'S in the limit a D — > 0. As a result, 9 a (x) ~ (2B'S)ln(\x — 
x~q\) in two dimension, and the corresponding 'electric field' energy cost ~ / d 2 x(V9 a ) 2 
diverges logarithmically as the size of system goes to infinity. 

The divergence of magnetic energy to order 0(B 2 ) indicates that the (classical) response 
of a non-magnetic impurity to external magnetic field is intrinsically non-linear and 
suggests that quantum effect may play an important role in determining the correct response 
of non-magnetic impurities to external magnetic field. In the following we shall derive in 
the continuum limit, the 1/S (spinwave) Lagrangian in the presence of background 9 S and 
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6 a fields, and shall show that the infra-red divergence in classical energy can be cured by 
formation of local magnetic moment around non-magnetic impurity once quantum effects 
are considered. 

To derive the spinwave Hamiltonian in the presence of magnetic field, we rotate our co- 
ordinate system locally on each site such that the local z— axis is always along the 'classical' 
spin direction. In this co-ordinate system, the Hamiltonian becomes 



j 



<i^io,j> 

—gfisB (SjSinOj + S\ x ^cos9i) — g^sB ^2(cos9jSj X ^ — sinOjS^). (5) 

The Hamiltonian can be rewritten in Schwinger-boson representation of spins in the usual 
way. In the large- S limit, we may write 

z vi = z t\ = z fi = z fi = \T( 2 S)i 

where Z(Z)fJs are spin a Schwinger-boson creation (annihilation) operators for site i on 
sublattice a. Expanding the Hamiltonian to order O(S), we obtain to order 0(B 2 ), 

H = E d + Hi/s, 

where 

<i^ioJ> L 

+ gfJ . B B £ {B^Zpti + g^ B B Ytfi)Zpfv ( 6 ) 

To further analyse our system we again go to the continuum limit. Introducing symmetric 
and anti-symmetric boson fields 

m = ^(zf(x)-zf(x)i 

7r(2) = -^(Zf(2)+Zf(£)), (7) 
and integrating out the ir(x) field, we obtain in the continuum limit an effective Lagrangain 



for <f>{x) field [ 10|,[lT| . Details of the technique can be found in ref. |10] and we shall not 



repeat them here. For x away from the impurity site (a < \x — xq\), we obtain in imaginary 
time 



L l/S ~ ]*r]*x± 



— -w a)V K — + e6 a )(J) + m 2 <J) + <J) + c 2 \V<J)\ 2 



(8) 



MS 

where e = B' and m 2 = S 2 B f2 /2. c 2 ~ 4d(JS) 2 is the spinwave velocity. We have set 
9 s (x) = B'S/2d and have neglected (V6 a ) 2 terms in deriving Lug- The later is of higher 
order (0(1/ S)) compared with the corresponding term in E c \. The most striking feature of 
the effective Lagrangian is that the 9 a field now appears as the r-component of a U(l) gauge 
field coupling to a charge boson field <ft, with the dynamics of the 9 a field governed by E c \. In 
particular, in the presence of B field, the non-magnetic impurity appears as electric charge 
generating a static electric field coupling to the charge bosons. Notice that the effective 
charge of the non-magnetic impurity changes sign when it is moved from one sublattice to 
another, indicating that the 'sign' of the charge is in fact a sublattice index |jlO| ,p 



The logarithmic divergence in 'electric field' energy in presence of non-magnetic impurity 
in E c i in two dimension will be removed if bosons of opposite 'electric' charge are nucleated 
from vacuum to screen the effective electric field, forming effective local magnetic moments 
around the impurity. The number of bosons nucleated from vacuum is ~ 2S, as can be 
seen easily by counting the number of 'charges' carried by the non-magnetic impurity. The 
magnitude of magnetic moment formed around the impurity is thus ~ S. The energy 
'cost' for nucleating the bosons can be computed by solving the corresponding Schrddinger 
equation for charge bosons moving in scalar potential of external charge of magnitude 2Se 
|T^| . To capture qualitatively the physics, we estimate the energy 'cost' by using a variational 
wavefunction ip(r) ~ Ae"'^'^ 1 , where £1 is determined variationally. Minimizing the energy, 
we find that ^ is of the order £ x ~ (JSc 2 /(g/is-B) 3 ) 1 / 2 ). Notice that £1 — ► oo as B — > or 
S — > oo, indicating that the formation of magnetic moment around non- magnetic impurity 
is a non-perturbative quantum effect which cannot be captured by usual spinwave theory. 
The energy 'cost' for nucleating the bosons is of order E\ ~ (2S)(m + 2SB' 2 ln(^ 1 /a )). 

Next we discuss the situation when there is finite concentration of impurities n randomly 
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distributed in the system, which is the case of experimental interest. In our continuum 
theory, the presence of finite concentration of random impurities is equivalent to putting 
finite concentration of 'charges' with magnitude 2SB' and random sign in our effective 
charge boson system. The behaviour of the system is very different from the one impurity 
case, since now electric fields originating from opposite charges will cancel, with remaining 
'electric field' energy of order ~ nE e , where E e ~ (2SB') 2 ln(l / a ) on average, I ~ rT x l d is 
the average distance between non-magnetic impurities. For E e > E\ ~ 2Sm or equivalently 
/ >> £ ~ a e Js ^ 29 ^ BB \ it is energetically favourable to nucleate charge bosons from vacuum 
to screen the electric field, and magnetic moments will be formed around non-magnetic 
impurities. The magnetic moments couple to each other weakly with J e jf ~ Je^ 1 ^ 1 implying 
that Curie behaviour will be present in magnetic susceptibility down to very low temperature 
~ J e ff- As concentration of non- magnetic impurities increases, / decreases until I ~ £ D . At 
this point a transition occurs where it becomes energetically unfavourable to nucleate bosons 
from vacuum to screen the electric field, i.e. there will be no magnetic moments forming 
around non-magnetic impurities at zero temperature when / < £ D . The boson-non-magnetic 
impurity bound states become excited states of the spin system! 

The boson-non- magnetic impurity bound states can still be observed as effective local 
magnetic moments at finite temperature T > m ~ g/iB-B if E e + 2Sm > E\ (corresponding 
roughly to ^ < I). At this energy range, the boson-impurity bound state appears as stable 
excited state of the spin system with excitation energy ~ (E 1 — E e )/2S < m. At finite tem- 
perature T > m, these states will be occupied by thermally excited bosons, forming effective 
free magnetic moments around non-magnetic impurities. It is important to emphasize that 
the formation of excited boson bound state around non-magnetic impurity is only possible 
when there is a gap m in the boson excitation spectrum. In the absence of non-magnetic im- 
purity, a direct computation of the spinwave spectrum indicates that the spinwave spectrum 
splits into two branches at small momentum q, with one branch carrying a gap larger than m 
in our effective Lagrangian, and the other branch gapless as q — > 0. The contribution from 
the splitting of the two branches to Hi/ S cancels to order 0(B 2 ), leading to our effective 
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Lagrangian with only massive bosons. The splitting of the boson spectrum reappears if we 
include terms to order 0(B 4 ) in H 1 / s , and the boson bound state around non-magnetic 
impurity will cease to be good eigenstate of the effective Lagrangian (the bound state turn 
into resonant state). However, the decay rate of the resonant state ~ B A is small compared 
with the energy of the resonant state ~ B 2 in the limit of weak magnetic field, indicating 
that the physical picture of boson bound state around non-magnetic impurity is still a good 
approximation to the system in weak field limit. 

As concentration of impurities increases and distances between impurities decreases fur- 
ther, the overlap between boson bound state wavefunctions at different impurity sites in- 
creases and a boson impurity band is formed (recall that the size of the boson bound 
state wavefunction is of order ^i). For I < £ 1; the bandwidth W of this impurity band 
is of order ~ (-y/m 2 + (c//) 2 — m) and the effective impurity potential strength is of order 
v ~ 2SB' 2 /(1/7) 2 . We find that the screening of effective impurity potential by charged 
bosons is weak and can be neglected. Estimating the scattering life time r we find that 
bosons with energy E < Ei ~ c//£ 2 are strongly localized in this impurity band using the 
criteria Er < 1. Notice that the number of localized boson states ~ rT l and decreases 
as n increases in this regime. The magnetic susceptibility will show Curie behaviour at 
temperature above TJ ~ Ei/ks- For T < TJ, the number of effective magnetic moments con- 
tributing to magnetic susceptibility decreases as temperature lowers. Roughly speaking, the 
contribution from non-magnetic impurities to the magnetic susceptibility can be described 
by a magnetic field and temperature dependent density of local moments n(B,T) which has 
the following properties: (i) n(B,T) ~ at the region qhbB » ksT, (ii) n(B,T) ~ n for 
I > £1, and ~ for I < £1 at the region k B T » gfisB + Ei, and is smoothly interpo- 

lating between the two regions. As temperature T — > 0,n(B,T) — > and the contribution 
to magnetic susceptibility from impurities is of order 5x ~ AS 2 nln{l/a ). 

Our theory can be extended to three dimension in a straightforward way. We find that in 
the weak-field limit, no magnetic moment is formed even in the limit of one single impurity 
at zero temperature because of absence of infra-red divergence associated with Coulomb 
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potential in 3-D. However, bound states of bosons can be formed at T > g^sB as in two 
dimensional case, and can still result in Curie-type behaviour in magnetic susceptibility for 
/ >> £1 ~ (c/gfisB) 3 . As concentration of impurities increases further, the boson bound 
states turn into impurity band and the effective number of 'free' magnetic moments n(B, T) 
decreases when temperature is lower than the boson impurity band bandwidth as in two 
dimension. 

Summarizing, in this paper we have carried out an analysis of the effects of replacing 
spins by non-magnetic impurities in Heisenberg antiferromagnets in magnetic field in two- 
and three- dimensions in a semi-classical 1/S expansion. We find within a continuum approx- 
imation that magnetic moments will form around non-magnetic impurities once magnetic 
field is put on the system, resulting in Curie-type behaviour in magnetic susceptibility when 
concentration of impurities is not too high, and temperature is not too low. Notice that our 
theory predicts that formation of local magnetic moments in presence of external magnetic 
field and non-magnetic impurities is a general behaviour of ordered quantum antiferromag- 
nets and is not restricted to particular materials. The theory explains the observation of 
Curie-type behaviour in magnetic susceptibility in La 2 Cui_ x Zn x 4: and Sr(Cui^ x Zn x ) 2 3 
compounds at temperature well below Neel temperature. In particular, the vanishing of 
local moments at T < gfisB is a theoretical prediction which can be tested experimentally. 
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